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A system of Brownian motions in one-dimension all started from the origin and conditioned never 
to collide with each other in a given finite time-interval (0, T] is studied. The spatial distribution 
of such vicious walkers can be described by using the repulsive eigenvalue-statistics of random 
Hermitian matrices and it was shown that the present vicious walker model exhibits a transition 
from the Gaussian unitary ensemble (GUE) statistics to the Gaussian orthogonal ensemble (GOE) 
statistics as the time t is going on from to T. In the present paper, we characterize this GUE- 
to-GOE transition by presenting the graphical expansion formula for the moments of positions of 
vicious walkers. In the GUE limit t — > 0, only the ribbon graphs contribute and the problem is 
reduced to the classification of orientable surfaces by genus. Following the time evolution of the 
vicious walkers, however, the graphs with twisted ribbons, called Mobius graphs, increase their 
contribution to our expansion formula, and we have to deal with the topology of non-orientable 
surfaces. Application of the recent exact result of dynamical correlation functions yields closed 
expressions for the coefficients in the Mobius expansion using the Stirling numbers of the first kind. 

PACS numbers: 05.40.-a, 02.50.Ey, 02.10.Ox 



I. INTRODUCTION 

Statistics of a set of one-dimensional random walks 
conditioned never to collide in a given time interval, say 
T, has its own importance in statistical physics, since, 
if we set T ^ oo, it realizes the one-dimensional Fermi 
statistics 0], and with T < oo it is used to analyze the 
models for wetting and melting phenomena We will 
refer to such non-colliding random walks and the contin- 
uum counterpart, non-colliding Brownian motions, sim- 
ply as vicious walks following the terminology used by 
M. Fisher . For thepioneering work on vicious walker 
models, see 0, 0, IE Hi 0- A generic setting of vicious 
walk problems is discussed in Q. Recently the interest 
on the vicious walks in mathematical physics is renewed 
and growing very rapidly, for close relationships of the vi- 
cious walk problem with the study of ensembles of Young 
tableaux and the symmetric functions [E 0|, the 
theories of orthogonal polynomials and random matrices 
^E and some topics of representation theory and 
probability theory 0, 0, 0| have been clarified. 

In an earlier paper Il7j . the continuum limit of non- 
colliding random walks on a lattice was taken by letting 
the temporal and spatial units At, Ax go to zero with 
the relation At cx (Ax)^ and a system of non-colliding 
Brownian motions was derived. Since each random walk 
tends to be a Brownian motion in this diffusion scaling 
limit, such a construction of non-colliding Brownian mo- 
tion is plausible, and indeed mathematical rigor can be 
established as a functional central limit theorem of vi- 
cious walks • The important fact is that the repulsive 



interaction among the obtained Brownian particles is no 
longer contact interaction as in the original vicious ran- 
dom walks on a lattice but is long ranged. The origin of 
this long-ranged interaction is the restriction of allowed 
configurations by the non-colliding condition, so that we 
can say it is an example of entropy-origin effective force. 
Moreover, the following setting was considered in our vi- 
cious walks; (i) even after taking the continuum limit, 
still we let the time interval, in which the non-colliding 
condition is imposed, be finite T < oo, and (ii) all the 
Brownian particles are assumed to start from an origin 
and the time interval of non-colliding is (0,T]. 

In this setting the process is temporally inhomoge- 
neous. At the very early stage t <^ T the repulsive in- 
teraction should be strong, since the Brownian motions 
will be restricted so that they will not collide for a long 
time period up to time T in the future. As the time t 
going on, the strength of repulsion is decreasing as is the 
remaining time until T, and attains its minimum at the 
final time t = T, at which there is no more restriction 
of motion in the future t > T. It was Dyson's idea that 
such systems of Brownian motions with long-ranged re- 
pulsion could have the equilibrium states, which can be 
described by the distribution functions of eigenvalues of 
random matrices in the ensembles appropriately specified 
by symmetry of the system |l8l| . The previous paper [T^ 
showed that the spatial distribution of vicious walkers 
at t ^ T realizes the eigenvalue-statistics of the Gaus- 
sian unitary ensemble (GUE), one at t — T does that of 
the Gaussian orthogonal ensemble (GOE), and the GUE 
to GOE transition is observed in the intermediate time, 
which is equivalent with the transition studied in the two- 
matrix model by Pandey and Mehta [T^ |20| . 
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In the present paper, we will characterize this transi- 
tion in distribution by calculating the moments of the 



2 



positions of N particles as functions of time t, 



MN.T{t,k) = 




k ^ 1,2,3,- 



(1) 



where Xj denotes the position of j-th vicious walker and 
( • )i the average at time t € (0, T]. In addition to such 
an interest of statistical physics, we can put emphasis 
that our present study possesses another importance as 
an interesting application of the graphical expansion the- 
ory. In high energy physics, graphical expansions for the 
matrix models of SU(A^) gauge theory were studied using 
ribbon graphs representing propagators and it was shown 
that the dominant graphs for large N are the planar ones 
and the leading corrections come from the graphs embed- 
ded in a torus, which are depressed by a factor with 
respect to the planar graphs 0,1121 . On the other hand, 
it was clarified that in the graphical expansion of SO{N) 
gauge theory the leading corrections only are depressed 
by a factor for large iV with respect to the domi- 
nant planar graphs, in which propagators can be repre- 
sented by twisted ribbons ^2^ |24| . In the gauge theory, 
the presence of non-Gaussian interaction terms of cubic 
or higher power is crucial and it leads to triangulations 
of random surfaces. Although the Gaussian matrix mod- 
els associated with the present vicious walker model arc 
not related with the triangulation problem of random 
surfaces and only give purely enumerative problems of 
surfaces, the proper structures of large N expansions in 
SU(A^) and SO(iV) gauge theories are also found in the 
GUE and GOE models, respectively. We can show that, 
if T ^ oo in our vicious walker model, the moments of 
the walker positions iQJ can be calculated by the graph- 
ical expansion method of GUE using the ribbon graphs 
and the results are given in the form of power series in 
the inverse of the square of matrix size N , 



ff>0 



iV2 



Here the coefficients eg{k) are the numbers of orientable 
surfaces of genus g made from 2fc-gon by some specified 
procedure [l^ (see also j2^ and Sec. 5. 5 in ^^). On the 
other hand, for the GOE we have to take into account the 
non- orientable surfaces as well as the orientable ones and 
the expansion is in the form of power series in \/N . In 
other words, in order to generate necessary surfaces, we 
need to use twisted ribbons, whose type of graphs is now 
called of Mobius graphs . For the moments (Q) of 

our vicious walkers, we will perform the Mobius graph 
expansion in the present paper. Now the weights of the 
graphs are depending on the time t; in the limit i ^ all 
the weights on the twisted ribbons are zero, but they are 
growing as time t going on, and dXt — T twisted ribbons 
are equally weighted as untwisted ones. This gives an- 
other characterization of the temporally inhomogeneity 
of the process and the GUE-to-GOE transition. 



Quite recently Nagao, Tanemura and one of the present 
authors applied the method of skew orthogonal poly- 
nomials and quaternion determinants developed for the 
multimatrix models to the vicious walk prob- 

lem and derived the quaternion determinantal expres- 
sions for dynamical correlation functions 1^3 • Using this 
result, we will present an expression of the coefficients in 
our expansion of moments using the Stirling number of 
the first kind. 

The paper is organized as follows. In Sec. II we briefly 
review the previous results reported in [T^ and give the 
precise definition of the moments which we will study. 
Graphical representations are demonstrated in Sec. III. 
Application of the result of [s^ is given in Sec. IV to give 
the expression for the coefficients of expansion and some 
concluding remarks are given in Sec.V. Appendices A and 
B are prepared to derive the expression of the density 
function used in Sec. IV and the 1/iV-expansions of the 
one-point Green function discussed in Sec.V, respectively. 



II. VICIOUS WALKS AND MATRIX MODEL 

We study a continuum model of vicious walks, the 
non-colliding Brownian motions in the time interval 
(0, r], constructed in [l^ [TtI I as the diffusion scaling 
limit of vicious random walks on a lattice. First we 
briefly review our previous results. The configuration 
space of the present N vicious walkers is R;^ = {x = 
{xi,X2, ■ ■ ■ , xn) S R^; xi < a;2 < • • • < xn}, where R 
is a set of all real numbers. The probability density of 
vicious walkers at time t S (0, T] with the initial condi- 
tion that all walkers start from the origin is denoted 
by pN^xit,^)- It was given as 



Pn 



T{t,^)^Ce-\^\"'^'hN{^)NN{T~t,^) (2) 



with c - 2-^/2r^(^-i)/4t-JvV2/j|^^^^r(j72) for X e 

where r(z) is the Gamma function, ft.jv(x) = 



R 



Y{i<j<i<N{^t^ Xj) and 



AAAr(s,x) 



dy det 

j^jv l<j4<N 



-(a:j-yf)/2s 



By usin g th e Harish-Chandra/Itzykson-Zuber integral 
formula |33l Is^ 35], we will see that pN,T(t,yi) is pro- 
portional to the integral 



(x)^ y dlJ j dA exp 



with 



■H{H, A) 



T 



2t{T-t) 



T 



t{T^t) 



HA 



rp2 



2t^{T-t) 



(3) 



A' 



where X is the N x N diagonal matrix with Xj£ — xjSji, 
and the integrals / dU and J dA are taken over the 
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groups oi N X N unitary matrices {U} and real sym- 
metric matrices {A}, respectively. The proportional co- 
efficient is determined so that the probability density is 
normalized. On the other hand, the integral over A in 
(O can be regarded as the convolution of the Gaussian 
distribution of complex Hermitian matrices H with vari- 
ance t{l — t/T) and that of real symmetric matrices A 
with variance /T , and thus we have the expression 

PN.rit, x) cx hN{^f I dU fiN.Tit, U^XU) (4) 



with 




(5) 



where and in the following we use the abbreviations 
and for the real and the imaginary part of the complex 



number z, respectively, i.e. z — z 
i = \/— 1. Remark that, if we set 



R 



t(2T - t) 



T 



for z E C with 



(6) 



N 

/ dUj2^ft^N,T{t,U''XU). 



(8) 



Next we introduce the integration measure for the N x N 
complex Hermitian matrices 

dH= H dHf, Yl dH],. 

i<j<e<N i<j<e<N 

Since dH (x dU x /iAr(x)^dx, if x = {xi, ■ ■ ■ , xn) are 
the eigenvalues of H and dx = O^Li dxj (e.g. see 20f ), 

and ^2^=1 ^f' ~ trH^'^ for H = XU with any unitary 
matrix [/, © with ||SJ) becomes 



(9) 



where ( / ) = j dH / fiN,T{t,H) for functions / of the 
elements of H with 



/2TTt 



n 



i<i<^<JV 



7rt 



e-iH],f/til~t/T) 



X TT , 



(10) 



and ~ 1 — t/T, c'^ ^N,T{t, cH) is equal to the probabil- 
ity density of the two-matrix model of Pandey and Mehta 
with the parameter a 0, |23| . Corresponding to chang- 
ing the parameter a from 1 to in the Pandey-Mehta 
two-matrix model, a GUE-to-GOE transition occurs in 
the time development of particle distribution in our vi- 
cious walks. 

Now we define the quantity, which we will study in the 
present paper; the moment of particle positions in the 
vicious walks. Since the distribution (|2J) of x is symmetric 
about the origin 0, all of the odd moments vanish. The 
even moments are defined and denoted as follows, 



Mjv,T(i, k) 




/ dx Vxf /9Ar,T(i,x) 
■'^< .7 = 1 



(7) 



for fc = 1,2,- • •. 



Note that 

where the sum is taken over all N'^'^ combinations of in- 
dices ji , j2 , ■ • • J j2k , and that Hji — i/j^ + iHjg with the 
Hermitian condition Hf- = Hjg,Hjj = the inte- 

grand ivH'^^ in ® is a polynomial of the N"^ independent 
random variables {Hf^; I < j <i <N}\J {H]f 1 < j < 
I < N}. Since the probability density H10|l is a product of 
independent Gaussian integration-kernels, we can apply 
the Wick formula with the variances 



{{H],f) = '^^'^'^\ l-5,,l (11) 

for 1 < j < i < N , where 5ji is Kronecker's delta. 
We can readily prove that (|11|) is equivalent with 



III. GRAPHICAL EXPANSIONS 

A. Wick formula 

First we notice that is invariant under any permu- 
tation of zi, a;2, ■ • ■ , xjv- Then |(7J| is written as 



MN,T{t,k) oc — / dx /lAr(x)" 



{Hjf Hmn) = {SjnSem + J SjmSen) , (12) 



where c is given by |(BJ and 



7 ^ 



2T -t 

The Wick formula for @ is thus 



(13) 
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il j2:---,i2fc ir6S2fc:R 

(-^jir(2fc-l) JT(2fc-l) + l '^J7r(2fc)i7r(2fc) + 1 ) ' 



(14) 



with the identification j2fe+i = ii, where the first sum is 
taken over all N'^'^ combinations of indices :/i , j2 , • • • , j2k , 
and the second one over the set of permutations S2k of 
{1, 2, • • • , 2k} with the restriction 

R: 7r(l) < 7r(3) < • •• < 7r(2fc- 1), 
7r(2j- 1) < 42j), l<3<k. 

The total number of the terms in the second summation 
is (2fc- 1)!!. 



B. An example: the fourth moment 

In this section, by performing calculation of the fourth 
moment MN^T{t,2) = (tri/^), we will demonstrate how 
to obtain graphical expansions from the Wick formula 



J 



p4(l with the variance H12|l. We start from the Wick 
formula for MM,T{i, 2), 



N 



jl J2 J3 J4 

-^Jl J2 ^i3i4 ) {-^j2j3 Hjijl ) {-^jlj2 Hjijl ) {-^j2j3 ^jaji ) | ) 



which has (2 x 2 — 1)!! = 3 terms in the summand of 
ji, ■ ■ • , j4. Substitution of (|12|l and binomial expansion 
give the 3 x 2^ = 12 terms in the form, 



/ 2 \ 2 12 

\ / — 1 



with 



Li 

is 

L9 
ill 



■^J4J4I 



^h33^0232 ^j3jl^j 
Jl J2 J3 J4 

^iii2^j2j3^i3ji^iiiil ^ 

h j3 di 

^jlji^j2j3^j2jl^j3j4J 

h J3 J4 

^jlj3^hU^j2h^j3Hl ^ 

12,33,34. 

^3131^3234^3234^3333 ^ 

h ,33 ,34 

^3134^3231^3234^3333! ^ 



^2 — ^ ^3l33^32j2^3334^343l1' 



31,32,33,34 



31,32,33,34 



31,32,33,34 



31,32,33,34 



31,32,33,34 



31,32,33,34 



^3132^ 3233^ 3334^3431! ' 

Jl J2,j3,i4 

ifi = ^3134^3233^3234^3331^^ 
31,32,33,34 

is = ^3133^ 3234^ 3234^3331! ' 

jl J2,j3,j4 

ilO = ^3131^3234^3233^3334! ^ 

31 ,32 ,33 ,34 

il2 = X! ^3134^3231^3233^33341"^ ■ 

r 



31,32,33,34 



Graphically, we prepare a square with four vertices la- 
beled ji , i2 , ja, ji in a cyclic order for each term as shown 
in Fig.l and connect the vertices ja and jb by a line for 
each Kronecker's delta Sj^j^ . We then regard these lines 
connecting vertices as the hems of ribbons connecting 
the two edges of the square. For example, the two lines 
connecting ji <-!■ J3 and j2 <-> J2 in the term L2 are con- 
sidered as the two hems of a ribbon, say ri, connecting 
the edges jij2 and j^js of the square, while the lines 
ja ji and ji are as those of a ribbon, say r2, 

connecting j^ji and jiji . There are two ways to connect 



two distinct edges by a ribbon, by untwisting as ri and 
by twisting as r2 in the above example, respectively. For 
each twisted ribbon, we put a factor 7. 

Next we take the summation over ji , • ■ • , ji in each 
term. Again consider the term L2 for example. Un- 
der the restrictions on indices specified by the Kronecker 
deltas, ji = js — ji, only two indices, say ji and j2, can 
be chosen arbitrary from {1, 2, • • • , N}. Then the sum- 
mation over all possible choices of indices gives N'^ for 
this term. The contribution of L2 is thus N'^j. We list 
up the contributions of all terms in Fig.l, and the sum 



L5 : N 



U: Ny 



L, : Ny 




FIG. 1: Mobius graphs 



of them gives 

Mw,T(i,2)/(cV2)2 
= {N^ + N^-f X 2 + N-f'^) x2 + N + N-/x2 + N^-f^ 

= A^' |2 + (47 + 7^)1 + (1 + 27 + 27^)^1 . (15) 

This shows that the 12 terms are classified into 6 equiv- 
alence classes, {Li,Lg}, {L2, L3, Lio, Ln}, {L^^Lu}, 
{is}, {^67 -^7} and {is}) with respect to the contribu- 
tion to the moment, and Fig.l implies that all graphs for 
the terms in an equivalence class are topologically equiv- 
alent. 



C. General formula 

For the general 2/c-th moment, M]^T{t,k), k > 1, 
we have {2k — 1)!! Wick couplings in the formula (|14|l . 
each term of which is the k products of the variances 
{HjiHjnn}- By applying (|12|l and expanding in 7, we will 
have the K = (2fc- 1)!! x 2*= terms, Mjv,T(i, k)/{c^l2f = 
Sfci demonstrated in the above section, one-to- 

one correspondence is established between terms {Lg] 
and graphs each of which consists of a 2fc-gon with its 
edges connected by k ribbons to each other. For each 
graph corresponding to L^, let Lpi be the number of 
twisted ribbons in the k ribbons and Vt be the "free in- 
dices" remaining after the identification of indices under 
the Kronecker delta conditions. Then the contribution 
from the term Li is given by A^^*7'^*. As mentioned at 
the end of the previous section, we consider the equiva- 
lence classes of the terms having the same contribution 
to the 2fc-th moment, and let each equivalence class be 
represented by a graph F. We let T^(F) and ipiT) be the 




for the fourth moment. 



numbers of free indices and of twisted ribbons. Moreover, 
we denote the number of elements in the equivalence class 
F by |F| . In other words, |F| is the number of ways to gen- 
erate graphs, which are topologically equivalent with F, 
using a 2fc-gon and k ribbons by gluing edges of 2fc-gons 
by ribbons. Define Q{k) be the collection of all graphs 
{F} generated by the present procedure. Then we have 



MM^rit, k) = 




Each graph F having no twisted ribbons, ^p{T) = 0, 
defines a way of drawing a graph on an orientable surface 
S'r, called a map, and each map specifies the surface S'r 
on which the graph is drawn (see, for example, H^H^). 
In general, the specified orientable surface has "holes" or 
"handles" and their number is called the genus g{Sr)- 
The genus is related with V{T), E{T) = k (the number 
of distinct edges; the original 2k sides of polygon were 
glued together in pairs by ribbons) and F{r) = 1 (the 
number of faces) through the Euler characteristic, 

x{Sr) = V{T) -k + l = 2- 2g{Sr). (17) 

Then the contribution from all graphs having no twisted 
ribbons is expressed as 

= (y) ^'^"^E^.w(ii) ' (18) 

where l{cj} is the indicator; l{a)} = 1 if the condition 
Lu is satisfied and l{a)} = otherwise, and eg{k) = 
J2reg(k) l{¥'(r)=o,y(r)=fc+i-2ff}|F|. 
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In the similar way, the graphs F having twisted rib- 
bons, i^(r) > 1, are considered to define non-orientable 
surfaces ^r- The genus g for non-orientable surface may 
be defined by the Euler characteristics as [s^l 

x{Sr) = 2^g{Sr) 

instead of H17|l . Then all the contribution to the moment 
from such non-orientable surface graphs is 



2 y ^ 

ree(fc) 

y) ^'^'E^j H^.mik)!-. (19) 

^ g=l ^ ^ m=l 



where Sg^m{k) = J2reg{k) l{y(r)=m,\/(r)=/c+i-3}|r|- The 
moment (|16|) is then given by the summation 

MN.rit, k) = M% j.{t, k) + Ml, r^{t, k). (20) 

It should be noted that 7 defined by is a mono- 
tonically increasing function of t and changes its value 
from 7 = to 1 as the time passes from t = to T. The 
above formula H19|l shows the fact that the contribution 
from Mobius graphs with twisted ribbons is growing in 
time and aX t = T twisted ribbons contribute with the 
same weights as untwisted ribbons (the GOE case). 



IV. CALCULATION BY DENSITY FUNCTION 



Set 



/OAr.T(i,x;T,y) = 



C 



-|x|V2t 



(iV!)2" 



ft.Ar(x)sgn(/lAr(y)) 



X det , , 

i<j,fe<jv \ ^27r(T - t) 



for X, y e R^. It is easy to confirm that PN,T{t, x; T, y) is 
invariant under any permutation of xi, ■ ■ ■ , xn and that 
of 2/1, • • • , ?/Ar, and 10) is equal to TV! / dy pN,T{t, x; T, y), 
if X g R^. Then the density function at time t is defined 
as 

pit, x)=N In dx, I dy pN.rit. x; T, y), (21) 



and the 2fc-th moment is given by 

MN.T{t,k) = / x'^^p{t,x)dx. 



(22) 



Let Hj{x) be the j-th Hermitian polynomial, satisfy- 
ing the orthogonality J Hj{x)He{x)dx ~ hjSji with 



hj = 2^ As shown in Appendix A, the general for- 
mula for the dynamical correlation functions of vicious 
walks reported in [23, HSl gives the expression 



p{t,x) = -e 
c 



J=0 ^ 



1 



(Ar/2 -T)! E + 2j + i),7^^^+2,+i(x/c). (23) 



3=0 



Substituting (|23|l into (|22|l and replacing the integral 
variable x hy y = x/c give 



Af-l 



e ^ dy 



iN/2 + j)\ 
2^-^y^{N/2 - 1)! + 2j + 1)! 

X / y'^''HN-i{y)HN+2j+iiy)e~y' dy 



1 

^(N-l)\^' 



(2fc)! 



22fc+^0F(Ar _ 1)1 ^ j!(2fc - 2j)! 



H2k-2j{y){HN{y)Ve-y dy 
H2k-2j {y)HN-i{y)HN+i {y)e^y^dy 



V (^/2+j)! . 
22fe+7v^i0f(jV/2 - 1)! {N + 2j + 1)1 



E 

i=0 



{2k)\ 



e\{2k~2iy. 



X J H2k-2i{y)HN-i{y)HN+2]+i{y)e ^ dy. 

In the second equality, we used the Christoffel-Darboux 
formula (see page 193 in p^) 



N-l 

E ir{HAv)Y 



3=0 



1 



2^0F(iV- 1) 
and the relation 



j-^ [{HNiy)y - HN-iiy)HN+iiy)] 



(2fc)! 



3=0 



mk-2j: 



-H2k-2j{y)- 



Now we apply the integration formula for the triple of 
Hermitian polynomials (see page 290 in ^|) 
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72''^/7^, if j + ^ + m = 2s is even 



r , I {s ^ - l)\{s - m)\ 

/ H,{y)He{y)H^{y)e-y dy ^ i and s - j > 0, s - £ > 0, s - m > 0, 



0, otherwise. 
Then we arrive at M7v,t(^, k) — Af^ fc) + Mj^ rp[t, k), with 

c^X'' (2fc)! A /fc\ / TV 



M^,At,k) - (^2j ^ ^ iN/2y.iN + i-k + ,mk-j-£-iy.ik-j + £ + iy.ji^ ■ ^'^^ 



Here we do not repeat the explanation how to character- 
ize the quantity eg(fc) in lfTH|l by using the formula (|^ . 
and only mention that it is obtained as the solution of 
the recurrence relation 

{k + l)eg{k) = {4k-2)eg{k-l) 

+ l)(2fc- l)(2A:-3)£g_i(fc-2) (26) 

with the boundary conditions 



1, if 5 = 0, 
0, otherwise. 



See m, Sec.5.5 in |2^ and ^ for details. 

In order to express the expansion in 1 /N for 125(1 , here 
we introduce the number s(n, k) defined as the coeffi- 
cients of the expansion 

n 

x{x — 1) ■ ■ ■ {x — n + 1) = ^ s(n, i)x^ 

for n > 1. It is known that s{n,£) x (— 1)"~^ is the 
number of elements in the set Sn of all permutations 
of {1, 2, • • • , n}, which are products of i disjoint cycles. 
These numbers s(n,£) are called the Stirling numbers of 
the first kind For example, there are 3! = 6 distinct 
permutations of {1, 2, 3}. We denote a permutation 1 
a,2 ^ b,3 ^ c simply by [a b c] (a, b,c £ {1, 2, 3}, a ^ b ^ 
c). We regard [2 3 I] as a cycle 1 ^ 2 ^ 3 ^ 1, [2 1 3] as 
a product of two cycles I ^ 2 ^ 1 and 3^3, and the 
identity transformation [I 2 3] as that of three cycles 1 
1,2^2 and 3^3. In this example, we see s(3, 1) = 2 
(for there are two elements [2 3 1] and [3 1 2] with £ — 1 
in 53), s(3,3) = 1 ([12 3]) and s(3,2) = -3 (other three 
permutations). For convenience, we will assume here that 
s{n,£) = ii n < 0, or £ < or n < £. Then we have 
expression (|19|l from (|25|l with the coefficients 



Ak) 



{2k)\ 



2' " j2^g+^ t 



It is easy to confirm that s(n, 1) = {—l)^~^{n — 1)!, 
s{n,n — 1) = —n{n — l)/2 and s{n,n) ~ 1 for any 
n = 1, 2, 3, • • •, and numerical tables of s{n, £) are found 
in 01 . For example, if we set k = 2, 1271) gives 
£1,1(2) = 4s(l, l)s(2, 2) = 4, £1,2(2) = s{2, 2)s(l, 1) = 1, 
£2.1(2) = 6{s(l,l)s(l,l)-|-2s(l,l)s(2, l)/3} 2 and 
£2^2(2) = -2s(2,l)s(l,l) = 2, and having £o(2) = 
2, £1(2) = 1, the result (HSJ is again obtained through 
the general formula (|2nj) with (^Hl and l(Tn|) . 



V. CONCLUDING REMARKS 

In the present paper we performed the graphical ex- 
pansions for the moments of positions of vicious walkers. 
The obtained formula with (^HJ and can be re- 
garded as the power series of the number of walkers N. 
Here we consider the large N limit. Let M^^xit, k) be 
the dominant term of H20() in 1. Then 



MN,T{t, k) 



N'+'eoik), 



(28) 



that is, only the contribution from the genus-zero ori- 
entable surfaces will survives in the limit N 00. It is 
well known that £o(A:)'s are given by the Catalan numbers 
Cfc (see, for example, |4^ ) 



1 f2k 



(29) 



xs{m, £)s{j — m + l,k — g — £ + 2). 



(27) 



k + l\k J' 
and their generating function is 

a(C) = ^{l-2C-yT^} 

00 

= ^C'C-fc. (30) 



fe=i 
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Corresponding to H28|l . define the density function p{t, x) 
as 



wfiere 



MN,T{t,k) = J x^''p(t,x)dx. 

Multiplying the both sides by z'^ and taking the summa- 
tion over k from to cx), we will have 



1 



1 - Vl - 2c^Nz 



1 — zx'^ 



dx, 



where (jJEl, IHI) and were used. If we set z = 

l/2c^N, it becomes 



p{t,x) 



-dx = 1. 



2N - {x/cY 
This integral equation can be solved as pC 
1 



-y/2N - (x/c)2, if < V2Nc 



0, 



otherwise. 



In the large N limit, the density function will keep a 
semicircle shape independently of the time evolution, in 
which only the width of the semicircle depends on time 
and simply scaled by c [s^- In other words, Wigner's 
semicircle law is universal in A'^ — > cxd. 

The universal property of random matrix theory at the 
large N limit and its finite- TV corrections have been stud- 
ied by calculating the Green hmctions in the form of 1/iV 
expansions in the field theory _21,.22„ .23., .43j. In order to 
compare our present results with the previous ones, here 
we consider the one-point Green function defined by 



GN,T{t, z) 




(31) 



with = Xj/{VNc) for z G C. It is nothing but the 
generating function of the moments 

1 °° / 1 \ 

We define the coefficients M!p^\t,k) in the expansion of 
the moment as 

/ 2 \ ^ oo 

M^,T(t,fc)= y TV'^+i^— M(")(t,fc). (32) 
Then we have the expansion of the Green function 



as 



oo 



(33) 



4")(i,z) = ^> '(-3) Mp{t,k) 



1 



fc=0 



1 



2z2 



(34) 



Since we have obtained the closed expressions for any mo- 
ments, (Uni with and ((JSJ, the Green function (EH) 

is completely determined, and the coefficients G^*'' {t, z) 
in its 1 /A^-expansion (|33|l can be derived for any order n 
through the formulae (|18|) , H19|) with the known result of 
eg{k) 25] and (|27|l . For example, as shown in Appendix 
B, the present results give the following expressions for 
the first three terms; 



G^°\t,z) 
G^^\t,z) 



1 + «(C) 



2z "(^.).x 



and 



^Yz'dC, 
1 



l-a(C)2 1-70(0' 
C d'^ 1 Q d 1 



(35) 
(36) 



3z 1 - a(C)2 6z 9C 1 - aiCY 
«(C) fo ^ A 7a(C) 



1 - 

«(C) 



2zC 1 - a{Qf 



27A-I 

97 J 1 - 7a(C) 

7a(C) 



3^1- _ 1 I (37) 
'97 J1-7«(C) ^ ^ 



with C = 1/(22^). By using H3U|I . we can show that set- 
ting 7 = 1 reduces them to Equations (24), (26) and (28) 
with /3 = 1 (the GOE case) given by Itoi, respectively, 
who derived them by solving the loop equations |4^ . On 
the other hand, we can confirm that, if we set 7 = 0, 
G^"V, z/\/2)/v^ gives the /? = 2 (GUE) results of Itoi. 
Some details are given in Appendix B. Our result is new 
and general, and it will reproduce the previous results 
of 1 /A^-expansions for the GUE and GOE as the special 
cases with 7 = and 1, respectively. 

Our formula also gives a power series in 7 
with A^-dependent coefficients. Let M{k;'^,N) — 
MN^T{t,k)/{c^/2f. Then 

k 

M{k;j,N)=Y,-/"'hMN), 

where fk,m{N) are the polynomials of degree fc + 1 in 
N for m — and of degree fc in A^ for 1 < m < fc. 
Explicit expression of fk,m{N) is immediately obtained 
from H25|) . We remark that the coefficients of the highest 
order in 7, fk,k{N), is equal to the zonal polynomials 
Z(k){si, S2, • • • , Sfc), if we set all the variables si — S2 — 

•••=sfe=Ar m. 

As shown by Q with ||SJ), the present system of vi- 
cious walkers can be regarded as a Gaussian model as 
matrix model, and thus exactly solvable as demonstrated 
in this paper. We would like to state, however, that the 
probability density Q of the positions of walkers x is 
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not in the simple Gaussian form multiplied by /lAr(x), 
when < t < T, due to the factor 7Vjv(T - t,x). 
Combination of the Pfaffian representation of this factor 
given in and the facts that Pf(A) = (det^)^/^ 

for any even-dimensional antisymmetric matrix A and 
det^l = e'''''^"^, (0) is written for even N as 

PN.rit, x) cx ft-Ar(x) cxp (-y(x)) 

with the non-harmonic potential 

IxP 1 



2t 



tr lni^(T-t,x; 



where F{s, x) is the N x N antisymmetric matrix with 
the element Fjk{s,K) = (2/-y/7r)Erf((a;fc — Xj)/2^/s) with 
Erf(tt) — due~^ . Then the present results for < 
t < T are nontrivial. 

In summary, we demonstrated the GUE-to-GOE tran- 
sition in time for the vicious walk model by presenting the 
graphical expansion formula of the moments of walker's 
positions. The weights of contributing graphs are time- 
developing and our formula interpolates the genus ex- 
pansion of orientable graphs for GUE and that of non- 
orientable graphs for GOE by using a time-parameter 
7 = t/(2T — t). The formula provides a power series 
in the number of walkers N and it was shown that the 
exact expression of dynamical correlation functions re- 
cently reported by Nagao et al. [s^ is very useful in 
order to evaluate the coefficients in the series. By com- 
paring with the previous results of 1 /N expansion of the 
one-point Green function, we showed that our results are 
general and valid. Further applications of the quaternion 
determinantal expressions of dynamical correlations of 
vicious walkers to the graphical expansions of more gen- 
eral types of moments {e.g. correlators of moments at 
different times l45j) will be interesting future problems. 



APPENDIX A: DENSITY FUNCTION 

Let Hj{x) be the j-th Hermite polynomial defined in 
Sec. IV. We can read the density function H21() from [s^ 
m by setting M = 1 as 

N/2-l 



(Al) 



where 

RfXx) = 7^/2^«,,H,(x/c)7-^/2, 



(A2) 



dyReiy) 



re 



X / dzdz' 

-OQ<Z<z' <QO 

9 , /r\4f+l 



p{y,z) p{x,z) 
p{y,z') p{x,z') 



(A3) 



with 



and 



a2ij = (c/2)2^(52o-, 
a2i+ij = {c/2f'+\52t+ij-m2i-ij). 



pix,y) = — X 



/2nt 



VMT - 1) 



The function p{x, y) can be expanded using the Hermite 
polynomials as 



p{x,y) 

exp ■ 



1-7 



V2TTty/2Tr{T - t) 

-(x/cf ^y'-/2T °° j/2 

.Y^l^H,{x/c)H,{y/VT). 



Then, for ^ = 0, 1, 2, • • •, we will have 



^2i{x) = — e 



j>2e+i 



7 



{]-(2e+i)}/2 



-Pj2e+iHj{x/c), 



'^2e+i{x) 



7 



0-2^)/2 



E 

j>2l 



where /3j^'s satisfy the relation 



-l3j2iHj{x/c), 



(A4) 



Ea 



= Sjs, < s < j. 



Substituting l|J2|l . (|X3|) and ((X4|) into ((JT|) gives (j^ . 

It should be noted that, though the M = 1 case of 
l3^ is equivalent with Pandcy-Mehta's two-matrix model 
[l9| . the present expression (|23|l is more useful for the 
moment calculation as shown in Sec. IV. 



APPENDIX B: COEFFICIENTS IN 
1/iV-EXPANSION OF GREEN FUNCTION 

By definition of M^\t, fc) in (jSU, lO and (dJ give 



k 

m— 1 



m>l 



k + m 



(Bl) 



(B2) 
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1 



12A:!(fc- 2)! 



m>l 



2k 



-(3m- 1) 



2k 

k + m 



(B3) 



where we have used H27(l . H29(l . and the fact ei(fc) = 
(2fc)!/{12fc!(fc - 2)!}, which is obtained from ^ with 
(1^ . Through (jSO)) and (|S3Il, (|BT|| immediately gives 
(1231 ■ In order to derive §^ and from and (|B3)| . 
respectively, we can use the identity given as Equation 



(C.5) in m, 



for ,„>0. 



n + m 



and its derivatives with respect to C. It is easy to see 
that 



and Gi^\t, z) is zero at 7 = 0. Moreover, we have ob- 
tained the following expressions for < 1 — 7 ^ 1 



G^^\t,z) 
G^it.z) 



^17^~2 z-^-2 



2(z2 - 2)3/2 



(l-7) + 0((l~7r), 



2z2 + 3 - 2zV'i2~~2 z(z2 + 6) - (z2 + l)Vi2~2 



4(z2 - 2)5/2 



4(z2 -2)3 



(l-7) + 0((l-7n, 



and 



4(z2 - 2)5/2 



at 7 = 0. 
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